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Abstract 



cn 

' We study slow collective motion at finite thermal excitations on the basis of linear response 

. theory applied to the locally harmonic approximation. The transport coefficients for average 

\ motion, friction 7, inertia M and the local stiffness C are computed along a fission path of 

■ ^"^^Th within a quasi-static picture. The inverse relaxation time /3 = 7/M and the effective 

. damping rate rj = 7/(2a/M|C|) are found to increase with temperature, but do not change 

I much with the collective variable. The values found for rj and /3 as well as their behavior with 

' temperature are in accord with experimental findings. 

0^ ' 
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'. 1 Introduction 



The nature of nuclear dissipation is not well understood as yet. For isoscalar modes at finite 
excitations the best information available presently comes from fission experiments, when comparing 
the decay rate of fission with the ones for emission of light particles or gamma rays. Nowadays it 
seems not only possible to deduce numbers for the effective damping rate 77 [|T|, but to gain 
information about its temperature dependence as well [Q. The authors of find an 77 which 
increases markedly with T, at small to moderately large values of the temperature. 

Such behavior is hard to understand for macroscopic models for which friction is either treated 
on the basis of two body viscosity or on the wall formula @], In the first case friction should 
decrease with T as T~^, in the second it would practically stay constant. Both results represent two 
contrasting pictures of the nuclear dynamics. The first one assumes collision dominance, and is thus 
related in a sense to the nuclear compound model or its classical analog, the liquid drop model. The 
second picture would apply if, between two encounters or "collisions" with the wall, the nucleons 
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are allowed to move freely inside the nuclear medium, as is suggested by the models of independent 
particles. In principle, the other two transport coefficients appearing in 77, namely M, C must be 
expected to change with temperature, too. But for the macroscopic models just mentioned this 
dependence will be weak. Therefore, the observed increase of r] with T gives strong evidence for the 
necessity of having a microscopic theory. One possible formulation of the problem and its solution 
is based on a specific application of linear response theory. For a detailed description of the latter 
we may refer to - ||^ and . In the present paper this theory will be applied to continue previous 
microscopic computations of the coefficients for friction 7, inertia M and local stiffness C, such as 
the ones of [^] and ||^. Emphasis will be put on their dependence both on temperature T as well 
as on the collective coordinate Q. 



2 The linear response approach to collective motion 

In this section we only want to outline briefly some basic theoretical features to have the most 
relevant formulas ready for explaining our computational analysis; we will follow largely [|10|, 0] 
and §. 

Similarly to the deformed shell model, we assume to have at our disposal a Hamiltonian 
H{xi,pi,Q) which depends on deformation through the shape variable Q. For the sake of sim- 
plicity we just take one such degree of freedom, which later-on will be used to parameterize motion 
along some given flssion path. However, it does not suffice to restrict this H{xi,pi, Q) to the model 
of independent particles, for which one would have (for A nucleons) 

A ^ ^ ^ 

Hsra{Xi,Pi,Q) =^h{xuPi,Q) ) (2-1) 

1=1 

where h{xi,Pi,Q) stands for the dynamics of particle 1. First of all, the expectation value of 
Hsm{£i,Pi,Q) does not represent the system's total energy. Secondly, as it stands this Hamil- 
tonian does not account for the effects of collisions. Following |]1T| the first deficiency is easily cured 



by adding to the operator part HsmiQ) a c-number term to get a HrmfiQ) as the representative 
for the (renormalized) mean field. This c-number term can be chosen in such a way that the static 
expectation value of HrmfiQ) (i-e. the one calculated for a time independent Q) contains the re- 
normalization terms of the Strutinsky method. (We choose the notation of |^ which differs slightly 
from [^, but the connection between both is readily established). The second issue can be taken 
care of by adding the effects of collisions when treating dynamical forces. In principle one might 
write 

H{Xi,p,,Q) = HrrrM) +V}^}{X^,P^) (2.2) 

with the restriction of having the residual two body interaction V^^^}{xi,pi) appear only in dynamical 
quantities like the response functions. How that can be done in practice will be explained below. 
Here we would just like to mention that we want to assume this incoherent interaction V^^]{xi,pi) 
to be independent of the collective coordinate Q. 

With these precautions taken into account we may say that the average {H{xi,pi, Q)) represents 
the total energy Etot{t) of the system, and now even in a dynamical sense. However, in case that 
the nucleus is isolated this energy must be a constant of motion. Hence the equation of motion for 
the Q{t) can be constructed from energy conservation by applying Ehrenfest's theorem |T2|. This 
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is to say we may write: 

= jE,^t = g(^^^ll^). = Q{F{x,,P^,Q))t (2.3) 

The remaining task then is to express the average {F{xi,pi, Q))t as a functional of Q{t). It is here 
that we shall exploit the benefits of linear response theory. The relevant operator F{xi,Pi,Q) is 
seen to be given by the derivative of the mean field with respect to Q; it is of a purely one body 
nature as the residual interaction V^eK^i^Pi) drops out when calculating this derivative. 



2.1 Local linearization 

Under certain circumstances, the task of evaluating the functional form of {F{xi,pi,Q))t can be 
considerably simplified. This will be so whenever the relevant Q can be handled as being close to 
some fixed value Qq. Quite naturally this would be the case for harmonic vibrations about a stable 
potential minimum at Qq. Fortunately, such a situation may be given even in the more general 
case. In this paper we restrict ourselves to study average dynamics, which means to neglect any 
statistical fluctuations in the variable Q, in which case it suffices to require collective motion to be 
sufficiently slow. Then Q will stay in the neighborhood of some given Qq for a ("microscopically") 
large time interval St; if necessary one may interpret the Qq as the Q{t) at a given time to- 

For any Q close to Qq one may evaluate the intrinsic quantity {F{xi,pi, Q))t by effectively using 
the Hamiltonian 

H{x,,p,,Q) = H{x,,p,,Qq) + (g - QQ)Fix,,p,,QQ) + ^{Q - Qq^ {^[Qq))^^^^^^ , (2.4) 

instead of the original one. It is obtained by expanding the Hamiltonian to second order and by 
approximating the second order term in a kind of unperturbed limit. By this we mean to evaluate 
the expectation value on the very right by a static density operator = Pqs{H{QQ)} which is 
determined by the Hamiltonian H{xi,pi,QQ) taken at Qq. Effectively, the only coupling term left 
between collective and intrinsic motion is then given by the term of first order in Q — Qq. It is not 
difficult to grasp the concept behind such an approximation: In this way global motion is described 
within a locally harmonic approximation. Here it was developed for average motion; a discussion 
of the general case can be found in [Q. One last remark on notation: In the sequel the F{xi,pi, Q) 
will be denoted by F whenever it is to be taken at Qq. 



Applying this concept to derive the equation of motion, it can be shown (see e.g. [0 and 
that eq. (|2.3| ) leads to the following form: 



/oo 
X{s)q{t - s)ds = 0. (2.5) 
-oo 

Here q = Q — Qm measures the deviation of the actual Q from the center of the oscillator ap- 
proximating the true potential in the neighborhood of Qq. The x is the causal response function 
associated to the dynamics of the nuclear "property" [F). It is given by 

x{t-s) = Q{t-s)ltr{pqs{Qo.TQ)[F'{t).F\s)]) 

= 2iQ{t - s)x"{t - s) (2.6) 
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The expectation value appearing here is to be calculated like the one encountered in (|2.4|) , namely 
by the density operator pgs{H{QQ)}. The same Hamiltonian appearing there is used to specify the 
time evolution in F\t). The pgs is meant to represent a thermal equilibrium at Qq with excitation 
being parameterized either by temperature or by entropy. The quantity k summarizes contributions 
of static forces which appear in second order. Anticipating the change in entropy to be quadratic 
in q(t), one gets for the coupling constant k: 



(2.7) 



with x(0) being the static response (the Fourier transform of the time- dependent response function 
(|2.6|) taken at frequency uj = 0) and x"'^ being the adiabatic susceptibility (sometimes x(0) is 



referred to as isolated susceptibility). It was shown in |10] that for temperatures not smaller than 
1 to 1.5MeV the coupling constant may effectively be calculated from the free energy /, or its 
stiffness d'^f/dQl, rather, using the formula 



k- 



T=const 



+ X(0) 



(2.8) 



2.2 Transport coefficients from the collective response 

Transport coefficients parameterize properties of a system whose time development is described by 
differential equations. The easiest way of getting such an equation from the integral form (|2.5|) is 
to expand the factor q{t — s) under the integral to second order in s [jl2|. In this way the common 
equation of motion of the damped oscillator is obtained 



Mq{t) + -iq{t) + Cq{t) = 



and the transport coefficients attain the forms: 

M = M(0) = -i r 



X{s)s^ ds = ^ 



2 diu^ I ^=0 



7 = 7(0) 



X{s)s ds 



C = CiO) = -x(0) - k^ 



or 

_ 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



They follow after evaluating the moments in time of the response function by Fourier transforms. 
The last equality in eq. (|2.12|) is a consequence of the expression ( p.8|) for the coupling constant. For 
obvious reasons the notion "zero-frequency limit" has been coined to portray the coefficients (|2.10 
1). 



We may note that the procedure just applied is borrowed from models where the basic equations 
of motion are of integro-differential type, and where the reduction to differential form commonly 
is referred to as Markov approximation. It is only because of the self-consistency underlying the 
derivation sketched above that in ( p.5| ) the common inertial terms are missing. As typically they are 
of order zero in the coupling they would make up differential terms from the start. Unfortunately, 
for the present case the method behind ( |2.9| - 12.12| ) does not always lead to a decent approximation 
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to ( |2.5|) . Let us mention just two problems. The most stringent one can be seen in the fact that the 
expression for the inertia, which is nothing else but an extension of the cranking inertia to the case 



of damped motion at finite temperature [|I2|, may become negative. Secondly, as both (|2.5| ) as well 
as ( p.9|) are homogeneous equations, the transport coefficients obtained in this way are defined up 
to a common factor only. To obtain all three coefficients, namely inertia M, friction 7 and stiffness 
C, one needs additional information. 

Following general concepts, this information can be obtained from a response function which 
parameterizes collective motion locally, and which will thus be denoted by Xcoiii^^)- It can be 
derived by introducing a (hypothetical) external force fext{t)F and by evaluating how the deviation 
of < F from some properly chosen static value reacts to this external field in linear order: 



6 < F >uj= —Xcoii{^)fext{^)- As shown in ||T0[ and the Xcoui^) can be brought to the form 



Xcou[u:) = . . 2.13 

1 + kx{uj) 



which is known to be standard for the case of zero temperature (see e.g. |T3| and [0). For finite 
excitations one needs an additional condition, namely that the motion is ergodic, in the sense 
of having adiabatic and isolated susceptibility equal to each other: x(0) = x""^- 

The dissipative part of Xcoui^) represents the distribution of strength over various possible local 
modes, which exhibit themselves as individual peaks. The corresponding "dispersion relation" or 
secular equation 

1 + X{^) = (2.14) 

is easily recognized as the Fourier transform of (|2.5|) . From (|2.14|) one gets a series of complex 
frequencies u^, which actually come in pairs uj^ = ±8^ — iT^/2. Each pair may be associated to the 



solutions of the secular equation of a damped oscillator like ( |2.9| ). Any q(t) satisfying ( |2.5| ) fully or 
correctly will contain information about all the Ui,. However, restricting oneself to just one such pair 
the integral equation effectively are reduced to a differential form. Since we aim at describing slow 
collective motion, it is natural to take the lowest ones, say uf. In addition to the two parameters 
specifying them, we have as additional information the strength of the corresponding poles. This 
feature enables one to evaluate all three transport coefficients M, 7 and C. In practice this can be 
done by replacing 

{Xcoll{i^)y^S < F >^ = -fext{^) 

4 (2.15) 

{Xosc{i^))-'5<F>^ = k\-Mu^-^iu + C)5<F>^=-f,^t{u), 

which is to say by approximating the response associated of the low frequency mode by that of 
a damped oscillator with modified transport coefficients. In practice the latter can be found by 
fitting the dissipative part of the oscillator response to the peak in the original strength distribution 

X'colli.^)- 

The appearance of the factor k'^ has a simple mathematical reason. Evaluating the contribution 
to the collective response from the two poles chosen, the numerator of (|2.13| ) has to be calculated at 
the poles' frequencies, which in the Xcoiii!^) leads to the factor xil^t) = —l/k. Please notice that by 
construction the oscillator response introduced here is the one for the quantity < F >t, while the 
transport coefficients written in ( p.l5| are those for the " Q-mode" . Indeed, the form p.l3| ) requires 
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that (on average) motion in F is related to that in Q hy k < F >t= Q{t) — Qo- Therefore, the 
corresponding transport coefficients must be related to each other like = k'^T, where T stands 
for M, 7, C. 

As just indicated, the coefficients introduced in this way reflect the structure of the collective 
response ( p.l3| ), for instance by way of the self-consistency between collective and intrinsic motion 



[| . Therefore, we will at times call these coefficients the self-consistent ones. It is possible to relate 
them to the zero- frequency limit. Suppose the lowest pair of frequencies cof lies sufficiently close to 
CO = 0. Then the tof may be obtained by expanding in (|2.14| ) the x(^) to second order in to around 



a; = 0. The result is easily recognized as the secular equation to (p.9|) with the transport coefficients 
being given by (|2.10| - p.l2| ). In some conditions have been derived for the zero- frequency limit to 
apply, for the case where the collective response function Xcoui^) is well simulated by the oscillator 



response function Xosc{^), i-e. where the reduction (|2.15|) does not imply any further approximation. 



From the equations shown in |^ it is easy to see that the following relations hold true: 

CiO) 

x(o) 



C = C(0)fl + ^) (2.16) 



Once more, the transport coefficients of the left hand side are the ones obtained for the ffi and those 
of the right hand side represent the zero frequency limit of the Lorentzian defined by the M, 7, C. 
Please observe, that the form ( p.l3|) may be turned around to express x{^) by Xcoii{^^)- Replacing 



the latter by Xosc{^) the functional form of x{^) is seen to be of Lorentzian type as well. 

Before we continue we would like to demonstrate that this concentration in a low frequency 
peak can be considered realistic for many situations, provided the temperature is not too small. In 
Fig.|l] we show for temperatures T = 1 and 2 MeV the strength function Xcoiii^) the potential 
minimum and at the saddle point, which later on will be called point A and C, respectively. (Details 
of the computation will be explained in the next section). Whereas at A the strength distribution 
reflects the typical behavior of stable modes, the one at C corresponds to unstable motion. For stable 
motion peaks are seen at higher frequencies, but the latter get washed out with increasing excitation. 
These features have been recognized and discussed before in [^] for quadrupole vibrations of ^"^Pfo. 



The resulting concentration in a low frequency mode is even more clearly seen at the instability, 
where it already occurs at the smaller temperature of T = IMeV. We should like to stress that 
this somewhat peculiar behavior of having the main peak appear at very low frequencies is strongly 
related to the fact that the motion is over-damped. This can be made quantitative for the case of 
the oscillator response introduced in ( |2.15| ). It is easy to convince one-selves that the position Um 
of the maximum of Xosd^) is found at 



= ^ ( 2\/l - signCr^^ + r/^ + signC - 2r]^ ) (2.19) 



^Please recall that for undamped motion (2.13) is identical to the response function of RPA one gets for the 
separable interaction {k/2)FF. 
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where w"^ = \C\/M. In the hmit of oo ujm turns to |C|/7 and the value of the response function 
taken at is identical to Xosd^^m) = 

Finally we like to put Eq-s.(|l| to |2l8| ) on more general grounds. Before they have been 
derived on the basis of the Lorentzian fit to the collective response function ( p.l3 ). However, quite 
generally, one may start with the coefficients of the zero frequency limit obtained from the full, 
original intrinsic response function. These coefficients come into play expanding again to second 
order in u around u = 0, this time, however, not the x{^) itself, but the {xcou{^))~^ = k + (^(ci;))^-'^ 
instead. In this way approximate versions of the self-consistent transport coefficients are obtained. 
It so turns out that they obey the same relation to the zero-frequency limit as given by ( p. 16 



2.18| ). Later on we shall compare numerically the M, 7,C of (|2.16| - p.l8|) with the corresponding 
coefficients obtained by the fit indicated in ( p.l5|) . We may indicate already here that typically 
C(0)/x(0) is very small and decreases with temperature. This can be inferred from the fact that 
C(0) drops with increasing excitation whereas x(0) turns out quite insensitive to changes in T 
(c.f.0]). For realistic cases, say above T ^ IMeV, C(0)/x(0) is of the order of several %. For this 
reason the C and the 7 of ( p.l(j| ) and ( |2.17D get close to their zero- frequency value. It is only the 
inertia M of (|2.18|) which may differ considerably from M(0), as a consequence of the second term 
in the second bracket. It is interesting to note, that at T = the ratio C(0)/x(0) might become 
of order unity if the spin-orbit interaction would not play a crucial role. This fact may be inferred 
from the analysis of quadrupole vibrations presented in sect. 8. 5 of (from eq. (8.5. 13b) it can be 
deduced that in this case C(0)/x(0) becomes identical to unity). 



2.3 CoUisional damping of nucleonic motion 

Let us turn to the damping mechanism used in our theory. From the discussion of (|2.2| ) one may 
anticipate that finally it is the residual two-body interaction V}^}{xi^pi) which causes damping, 
first on the microscopic level, and then by way of self-consistency for the collective motion as well. 
Because of the assumption of this interaction being independent of Q it enters the game only through 
the Hamiltonian H{xi,Pi,QQ) which explicitly appears in the response function ( p. 61) . To evaluate 
this expressions fully for some given V}^} would be too a tremendous task. We therefore use a scheme 
which borrows from the way one would treat the effects of collisions in time dependent mean field 
theories like ETDHF or its classic versions as given by the BUU or Landau- Vlasov equation. In this 
paper we will just state the final expressions referring both to as well as to earlier publications 
(for detailed lists see and [||). 

The Fourier transform of the dissipative part of the intrinsic response function finally writes like 




Here, n{x) is the Fermi function determining the occupation of the single particle levels \k >. The 
latter are the eigenstates of the Hamiltonian /i(x, p, Qo) with corresponding energies e^. The gk{uj) 
represents the distribution of the single particle strength over more complicated states. It is here 
that the effects of collisions come into play. Neglecting them, which in our present language means 
to say putting V^}^} equal to zero, the ^^(a;) would simply be given by Qk{(^) = 27iS{huj — e^). 
Conversely, a finite V^^^] gives reason for finite self-energies for which both real and imaginary parts 
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are considered according to the formulae T,{u!,T) = T,'{uj,T) — iT{u,T) /2, with 

To 1 + ^ [{nuj - fiY + n^T^l 
and /i being the chemical potential. Then the gk{uj) becomes 

Qkiuj) = (2.22) 

In ( [^.21D the I/Fq represents the strength of the "collisions", viz the coupling to more complicated 
states. The cut-off parameter c allows one to account for the fact that the imaginary part of the 
self-energy does not increase indefinitely when the excitations get away from the Fermi surface. In 
the present calculation we choose Fq ^ = 0.03MeV~^ and c = 20MeV. We would like to draw the 
readers attention to the frequency dependence put into (|2.21|) . Speaking of the collision term in 



BUU type equations once more, the latter would have to account for subtle memory effects if it 
should cope with the form we use for the self energies. 

So far it was not specified whether or not the sum over j, k should include diagonal matrix 
elements. Quite generally, they measure quasi-static properties of the system and they are respon- 
sible for the "heat pole". The latter shows up at = either in the relaxation function or in 
the correlation function associated to the dissipative part of the response. It was argued in to 
neglect contributions from this heat pole, last not least to force ergodicity to be given. Within our 
model this implies to restrict the summation in eq.( |2.2(]| ) to non-diagonal matrix elements. Such a 
restriction has thus been done also in the present computations. A more elaborate discussion will 



be published in ||T6|. 



2.4 Friction in zero-frequency limit 

Because of its great importance we like to address specifically the calculation of the friction coeffi- 
cient, in particular its temperature dependence. The origin of the latter can be made transparent 
for the zero- frequency limit. As we shall see below, numerical evidence tells one that this limit 
represents the actual value quite well for not too small temperatures (c.f. also 0, [0] and |TB|). So 
let us evaluate 7(0) by inserting eq.( |2.20| ) into eq.( p.ll|) . One obtains: 



jk 



This expression can be simplified further using the Sommerfeld expansion for n{Q). To leading 
order one gets 

7(0)^|^El^.fcl'^?.(/^)^^(/^) (2-24) 

jk 



Here, the spectral density Qk{l^) is to be calculated at the frequency hoo = jj,. According to ( |2.21 
this implies to use for T{^,T) a constant value, independent of u. For such a situation it does not 
make much sense to use refinements like the cut-off in the frequency dependence. Moreover, as the 
real part of the self-energy vanishes at hu = fi (see 0) we may use 



r(/i, T) = ^n'^T'^ ^ 0.3T2 (2.25) 
Fq 



8 



and 

- (e.-,;^:r(..r) 

It may be noted that this approximation is very much related to the typical relaxation time ap- 
proximation (see 0) for the collision term of the BUU or Landau- Vlasov equation. 

Finally we may note that for the simple model outlined here the dependence of 7(0) on T 
appears to be synonymous to that on a constant damping width F. The validity of the Sommerfeld 
expansion has been studied numerically by computing eq.( |2.23 ) for several temperatures, and for 



such a model of a constant F instead of T{uj,T). It was found that the temperature dependence 
due to n{Q) is weak, which in a sense justifies the use of ( p.24| ). Therefore, a simple estimate of the 
temperature-dependence of friction may be obtained by using a constant width F and by relating 
the latter to T by way of ( p.25|) . The dependence of 7(0) on F may be studied explicitly within a 



simple model. Take a schematic nucleus consisting of A/2 protons and A/2 neutrons, which occupy 
oscillator shells with principal quantum numbers 1,2, ...,A^ and neglect the spin-orbit coupling. 
For quadrupole vibrations around a sphere one may proceed along the lines described in sect. 8. 5 of 
rT| . Applying this model to ( |2.24| ) with ( |2.26| ) one gets 



4T^{2huJo)'^(T 1 1 

97i[{hLJof + F2] + F2 + (5/icuo)2 + F2 + ■ ■ ^" 

The quantity a represents the dependence of the average squared matrix elements on A^. To leading 

'3 

-2^ 



order one gets a ~ (fA)3 (c.f.lll). 



2.5 The single-particle Hamiltonian 



For the present computation the single particle Hamiltonian h{x,p, Q) of ( p.l|) is chosen to be given 
by the two-center shell model: 



2m 



+ V{p, z) + Vis{x,l, s) + Vi2{x,l) 



(2.28) 



This model has been developed some time ago [|T^, (see also [|^) and has been applied successfully 
to heavy ion collisions and fission. It is comparatively simple numerically but nevertheless is able 



to describe quite well [|T8] the single particle energies and wave functions near the Fermi level. 

Here V is the two-center potential in cylindrical coordinates with m, s and / being the nucleons' 
mass and the operators for spin and angular momentum, respectively. A detailed description of the 
construction of V can be found in [|T7| and As basic shape parameters there are the distance 

= Z2 — Z\ between the centers Zi of the two potential wells, a neck-parameter e, a mass ratio 
a and deformations 5i, bi of the two fragments. To demonstrate the geometrical meaning of the 
shape variables, we show in Fig.^ the potential along the z-axis and the associated nuclear shape. 
As usual the nuclear surface is identified as an equipotential one of the two-center potential, and 
volume conservation is required for the uniform density inside the surface corresponding to the 
Fermi energy. The volume then has a value given by the sphere with radius i?o = tq^^'^^, with A 
being a mass number of the nucleus. In the present calculation, we take the radius parameter tq to 
be equal to 1.2 fm. The associated oscillator frequency Ulo^ for the spherical shape is then given by 
fiwo = 41A-V3. 
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The momentum-dependent part in eg. ( p. 281) consist of a spin orbit-coupling term Vis and an 
/^-term Vp. For them the angular momentum / is described in the stretched coordinates. The 
strengths Hi of Is— and KjyUj of P— are taken from |]18|. For a given shape the Hamiltonian (|2.28|) is 
diagonalized within the basis of eigen-functions constructed for the the two-center potential without 
a "neck-correction term"; hereto states up to energies of {NQ + 3/2)huo are taken into account, with 
A^o = 20. The "neck-correction term" is used to remove the cusp which naturally arises at z = 
if one just puts together two oscillator potentials (see Fig.|^; for more details we refer to [0 [p!9| ). 
The deformations 6i of the fragments are described by the ratio of the semi-axes and bi as 
ai/bi = + 2/36i) / \J{1 — 4/35i), with i = 1 and 2 for the left and right fragments. The neck- 
parameter e is defined as the ratio e = Vq/V with Vq and V being the heights of the potential 
barrier at the origin, calculated with and without the "neck-correction term", respectively. 

As mentioned earlier the transport coefficients shall be computed along a "fission path" only. 
In this paper we want to identify the latter with a line of minimal potential energy. This energy 
shall be approximated by the one of the cold liquid drop, calculated for a two dimensional subspace, 
which comes about in the following way: We concentrate on symmetric fission, for which the mass 
ratio a becomes unity, and choose the deformations of the left and right fragments to be equal to 
each other, which means to have = ^2 = 5. Finally, the neck parameter, which is sensitive to 
specifications of the scission configuration, is fixed to the value e = 0.4. In a previous study 
such a choice of the "fission path" has been seen to reproduce well the observed average kinetic 
energy for thermal fission of ^°°P6. The family of shapes from which the "fission path" is evaluated 
is defined by the two parameters zq and 5, and the motion along it shall be parameterized by a 
variable ri2. It is defined as the ratio ri2 = Ru/i^Ro), where Ru = R2 — Ri measures the relative 
distance between the centers of mass of the two fragments (see Fig.|]), and Rq stands for the radius 
of the spherical configuration. For this variable the single-particle operator F is defined as 



dh / dri2 95 dri2 
dri2 \ dzQ dzo d5 



-1 



' dh d6 dh ' 
dzo dzQ 86 



(2.29) 



where the derivative d6/dzo is to be taken along the fission path 6 = S{zq) 



2.6 The deformation energy 



For the collective response function eq.( |2.13D , we need the coupling constant k as given by eq. 
This expression can be evaluated from the free energy / and the static response x(0)- The free 
energy /(ri2,T) will be written as a sum of Coulomb and surface energies plus the shell correction 
part: 

/(ri2,T) = fcoui{ri2,T) + /.„./(ri2,T) + /.e(ri2,T). (2.30) 
The free Coulomb and surface energies is approximated by 



f Coui{r 12, T) = fcoui{ri2,T = 0)(1 - aT^ 



and 



(2.31) 



(2.32) 



f SUV fir 12, T) = fsurfiri2,T = 0)(1 - /?T^) 

with the values of a and /3 taken to be 0.000763 and 0.00553 MeV-^ 
fcouiiri2,T = 0) and fsurf{ri2,T 
there are several sets of parameters. Here, we choose the set with the values of a = 0.65/m, 



I|. The free energies 
0) have been evaluated according to [^. In Table III of |2^ , 
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tts = 21.836MeV", Kg = 3.48 corresponding to the radius parameter tq 
correction we assume the form (c.f.[|l^]): 



1.2 fm. For the shell 
(2.33) 



fsciri2,T) = fsciri2,T = 0)r/ sinhr 
with r = {2'it'^T) / (huo) . The shell correction fsc{ri2,T = 0) is evaluated according to [23 



3 Numerical results 

We want to evaluate transport coefficients along the fission path for symmetric fission of ^^^T/i and 
like to study their dependence on temperature T and their variation with the shape variable. As 
previously mentioned, the fission path will be defined through the static energy. This is justified in 
case that the collective motion is strongly damped, which implies that the system "creeps" down 
the potential landscape from saddle to scission. Indeed, later on we will find this hypothesis justified 
by the numerical values of the transport coefficients. 



3.1 The static energy and the stiffness coefficient 

In Fig.^the hquid-drop energy Vldm{zq, 6) = fcouiiT = 0) + fsurfiT = 0) is shown as function of zq 
and 6. The curve along which this energy is minimal for each zq is given by the long-dashed curve. 
It connects the point A of the spherical shape with the saddle point C and then follows down to the 
scission configuration. The latter part refiects the path of steepest decent. This trajectory is taken 
as our "fission path". It is clear that such a static path will change with temperature, because of 
the dependence on T of both the shell correction as well as of the macroscopic free energy. Such 
refinements are neglected here, for the following reasons. First of all, we concentrate on excitations 
where shell effects are largely gone. Secondly, we know from numerical experience that the path 
of minimal potential energy, as calculated from the "macroscopic part", is not very sensitive to 
changes in T. 

In Fig.§, we plot the free energy /(ri2,T) as a function of the shape variable ri2 for different 
values of T along the fission path shown in Fig.^. We restrict ourselves to T > IMeV and neglect 
pairing correlations, which should be taken into account for quantitative studies at lower tempera- 
tures. The second derivative of the free energy C(0) = d'^f/dQl is plotted in Fig.|^ as function of 
T. At not too small temperatures the actual stiffness C turns out to be close to this zero-frequency 
limit; ( |2.16| ) and the discussion given there. 



3.2 The friction coefficient 

a) A schematic study of the zero frequency limit 

We first would like to work out a few general features at the example of the zero-frequency limit 
(|2.24|) applied to quadrupole vibrations. In Fig.^ the result of a computation is shown for the doubly 
closed shell nuclei ^^^Yb and "^^^Pb, neglecting spin-orbit coupling in the first case. For the width 
the formula r(/i, T) = 0.3T^ is used and as single particle model a harmonic oscillator is taken with 
spheroidal deformation 6. For small deformations this 6 is related to 02 by 5 ~ (15/87r)a2, where 
a2 is the common parameter appearing in (i? = -Ro(l + tt2^2o(^) V'))) |13[, |Tl|]. 

Fig.p shows that at small T friction increases strongly with T, then it reaches some maximal 
value at some Tmax and decreases afterwards. Taken simply as function of the width F, such 
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behavior may be inferred from the semi-analytical formula ( p.27|) . With respect to the dependence 



on temperature one must be a little careful, as the result shown in the figure was obtained for 
r(yU, T) = 0.3T^, neglecting the influence of the cut-off parameter c in ( p.21| ). A finite cut-off 



parameter would weaken the T-dependence somewhat (see P), but it should be clear that the 
present model study only aims at a qualitative understanding. After all we are neglecting subtle, 
but important details of collisional damping like the width's dependence on frequency. The value of 
the width T^ax at which 7(0) reaches its maximal value is strongly related to the energy differences 
\ej — Cfcl of those single-particle pairs {j,k) which give significant values to the transition matrix 
element \Fjk\'^. From eq. (|2.27| ) one may deduce this Tmax to be of the order of Huq in case of no 



spin-orbit coupling. With spin-orbit couphng the values of Tmax will get smaller, since then there 
will be important contributions from pairs (j, k) which lie closer in energy. The corresponding T^ax 
may be estimated to 0.3 x T^^^ ~ Tmax- For example, it can be seen that the values of T^ax are 
6.7 Mev for ^'^^Yb and 2.5 MeV for ^ospfe. 

Let us address now the value of the friction coefficient and compare it with the wall formula. 
The maximal values of 7(0)/^ are seen to be 40 and 180 for ^'^^Yb and "^^^Pb, respectively. These 
numbers reflect a general behavior in that friction can be said to become larger if the spin-orbit 
coupling is taken into account, for given mass number. However, altogether the magnitude is much 
smaller than that of the wall formula, for which one would get 530 and 900 for the two cases 
considered, namely ^'^^Yb and "^^^Pb, respectively. The latter numbers can be estimated as follows. 
For the spheroidal deformation 6 the wall friction 

3 r 

Iw = -,pVF f vldS (3.1) 



can be shown to reduce to the simple expression 7^//i = 2! vOttAs ^ 0.76/1 3 (p is the particles 
density, vp - Fermi velocity, Vn - wall velocity normal to surface). 

The small value obtained for our friction coefficient is related to the "gap" in the single-particle 
spectrum, which in the expressions for the friction coefficient appears in the denominators (mind 



(|2.27|) ). For doubly closed nuclei the contribution to friction comes from states whose energy 
difference \ej — ek\ is rather large. This energy difference becomes especially large for the spherical 
oscillator potential (for which only states with \ej — ek\ = 2hujo contributes to response function 
and consequently to the friction). 

b) General results 

Let us turn now to discuss results which are obtained using our complete theoretical input, 
namely the realistic two center shell model, applied to expression ( |2.20| ), with the collisional damping 
as given by (|2.22|) and ( p.21|) , together with the realistic coupling constant. The coefficient 7 is 



shown in Fig.|^ as function of temperature, for the points A and C along the fission path indicated 
in Fig.^ For comparison on the right the values of the so called modified wall formula are shown 
as well, which reads '~im.w = 0.277^ 



Two computations are presented, one based on the fit of the oscillator response (see ( |2.15| ) 
and one based on the formula ( p.l7|) (with C(0),x(0) and 7(0) being calculated from the original 
x(ci;), without employing any fit). It is seen that both results agree almost perfectly well for all 
temperatures used at those points where the local stiffness is negative, and where, hence, the local 
motion is unstable. At the stable point A there is some disagreement at smaller temperatures, but 
even at T = 1 MeV the difference is smaller than about 30% and at T = 1.5 MeV it is almost 
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negligible. The origin of the deviation of these two computations from each other can easily be 
understood looking at the collective strength distributions discussed above with the help of Fig.[l|. 
It is only for small temperatures that this distribution cannot be approximated by one Lorentzian. 
We should like to note that the friction coefficient given by ( p.l7| ) is almost identical to the one of 
the zero-frequency limit. The correction given by (C(0)/x(0))^ only amounts to about 4%. 

Finally, we exhibit in Fig.|^ the coordinate dependence of the friction coefficient for various 
temperatures. It so turns out that for the present case the result at T = 1 MeV is very well 
represented by the modified wall formula. 



3.3 Mass parameter 

Finally we turn to the inertia for which some new problems appear, as may be demonstrated with the 
help of Fig.|^. There the temperature dependence is exhibited for three computations of the inertia, 
performed at the point A. Firstly, there is the one from the fit and secondly the one calculated 
from the expression ( |2.18| ). Like for the case of friction, both agree with each other for higher 



temperatures, say above 1.5 MeV, and for the same reason. Above this temperature essentially no 
higher modes exist, at least not for the present model. It is important to note that, different to the 
case of friction, for the inertia the zero-frequency limit M(0), as given by the "cranking inertia" 
(|2.10|) fails to be a reasonable approximation. Looking at (|2.18|) it is only by way of the second 



term in the second bracket that the M calculated this way is close to the one of the fit. As has 
been noted before, the zero-frequency mass may become very small, often even negative. 

Let us address now the third curve of Fig.^ which is marked by " sum rule" . As is well known 
from the case of zero damping, see e.g. or |jll|, there should exist the following, general relation 



of the energy weighted sum to the inertias of the various possible modes with frequencies cj^ (see 
also §): 

S={[F,[H,F]]) = - r dcoxMu; 

7T J-oo 

= Ev/ ^^^otn^ V = E TJT^ (3-2) 

(Here, in the last expression again the factor fc^ appears because we take the inertia of the Q-mode 
but relate it to the sum of the F-modes). As may be inferred easily from this formula, the inertia of 
one mode cannot be smaller than h'^/{Sk'^). Unfortunately, this is not necessarily the case without 
further precautions. 

Essentially two problems appear here. First of all, evaluating the integral 

S = 2— duxcoui^)^ (3-3) 

TT JO 

as function of the upper limit ujmax one realizes bad convergence. This fact is related to the 
construction of the response function, particularly that of x"- The forms given in ( |2.22D and ( |2.21| ) 
cannot be expected to be good at large frequencies. The second problem comes up whenever the 
single particle potential does not warrant the density of the shell model to agree sufficiently well 
with the one which is used in the definition of the shape variable. Usually, this is all right for the 
simple model of the deformed oscillator, for which it is easy to prove that the inertia corresponding 
to the total sum, i.e. h'^/{Sk'^), is given by the Mirrot irrotational fiow. Indeed, for this case it 
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has been possible in to show that the self-consistent inertia M(uji) turns into Mirrot for large 
temperatures. Again, the very fact that M{uJi) must approach the value given by the sum rule 
follows from the observation that with increasing excitation the strength concentrates in the low 
frequency mode. 

In a forthcoming paper [|l^ we will try to solve the second problem mentioned above, which 
involves the relation of the density distribution to the shape of the potential. For the present case 
we suggest the following pragmatic procedure. We just fix the sum rule value by the inertia Mirrot 



of irrotational fiow, which in turn is evaluated by the method of Werner- Wheeler To put it 
differently, we scale the inertia M = M{uJi) for the low frequency mode such that it becomes equal 
to Mirrot at large T. At smaller temperatures this M will attain larger values, of course. The result 
of such a computation is shown by that curve which in Fig.^is marked by "sum rule". 

In this context it is worth recalling that at somewhat larger temperatures collective motion gets 
over-damped such that inertia does not play a big role anymore. It is reassuring therefore to see 
that for such cases the fit of the oscillator response does not depend very much on the inertia (mind 
the discussion given above next to ( ^.19| )). 



In Fig. 10 we address the coordinate dependence, as obtained for various temperatures. From this 



figure one may observe that above T = 2 MeV the values do not change anymore with excitation. 
3.4 The inverse relaxation time (3 = ^/M 

In the literature the quantity P = ■y/M has come into use (see e.g.|[^ to [^) for parameterizing 
the strength of the friction force, where for M commonly the irrotational fiow value is taken. The 
inverse quantity (3~^ has a dimension of time and actually measures physically the (local) relaxation 
time to the Maxwell distribution in collective phase space. Taking our results for friction and inertia 



from Figs.p and [1^ we get a coordinate dependence of P as shown in Fig.^. It is seen that this 
dependence is much weaker than the one of 7 and M themselves. This is to be expected in a sense, 
as any geometrical factor contained in the individual quantities drops out when building the ratio. 
Actually, the variation along the fission path is weaker than the one with excitation. Incidentally, 
it is seen that /3 decreases somewhat with deformation. 

The values of (3 seen in Fig.|ll]are in the range of those published earlier for quadrupole vibrations 
of "^^^Pb [|15| 13], having the tendency of being slightly smaller than in this previous case, which is 



no surprise in the light of the discussion presented above for friction. In any case, these values are 
comparable to those associated with "linear response theory" in the compilation given in Fig. 13 of 
0. As seen from Fig. 12 of the same paper, they are in reassuring agreement with numbers found 
adequate in many theoretical descriptions of fission accompanied by emission of light particles. 

In Fig.0 our results for f3 are compared with that suggested by Frobrich and Gontchar for a 
phenomenological description of fission of the type just mentioned Their deformation parameter 



q is the same as our parameter ri2. It may be said that their model has also been quite successful in 
explaining observed experimental data, but quite apparently, their (3{q) is in violent contradiction 
to our picture. In case it will turn out that their model is safe otherwise, leaving no room for 
modifying the /5(g), one might be inclined to conclude that the friction necessary in this case must 
be of different physical origin. One possible candidate would be a dissipation mechanism which 
comes from the so called "heat pole" and which in [§] has been seen to be identical to the one of 
"diabatic" motion [E^. This would be an interesting perspective as it might indicate that motion 



beyond the saddle is "fast", even at high excitations, in contradiction to the general believe. 



14 



3.5 The effective damping rate 77 = 7/(2yM|C|) and the dynamics of 
fission 



For a damped oscillator the quantity rj = •y /2^M\C\ indicates whether motion is under-damped 
[t] < 1) or over-damped [r] > 1). Therefore, it is a good measure for the effective degree of damping, 
except around points of inflection where the local stiffness becomes zero. Furthermore, it is exactly 
this quantity, if calculated at the saddle point, which according to Kramers' famous formula 



Rk = (^^1 + Vs- Vs) ^ e-^/^ ^ [^^l + 77, - r^, j Rbw (3.4) 

determines the deviation of the fission decay rate Rk from the Rbw of the Bohr- Wheeler formula 
{zu^ = \C\/M). Within the linear response approach t] has been calculated as function of T before 
(see [^]), but only for the case of quadrupole vibrations around the spherical minimum of lead. 
From all the information presented above we are now in the position of calculating rj along the 
fission path. However, in the light of formula ( |3.4| ) it may suffice first to study it around the barrier. 
Indeed, this formula applies if, from the full dynamics across the barrier, it is only the motion 
around the stationary points which matters. The latter situation can be seen to be given if the 



height B of the barrier is much larger than temperature: B /T ^ 1 (|^, see also |29|). Effectively, 
one may then approximate the true barrier by a model potential which consists of one ordinary 
oscillator for the minimum plus an inverted one for the barrier, joined smoothly to each other (see 
p9|). Unfortunately, often ( p.4| ) is applied even in cases that this basic condition on T and B is 
violated. 

To study these questions for our example we have calculated all relevant quantities by performing 
averages both around the potential minimum as well as around the barrier. For the first case this 
means to average in the regime from = 0.375 to some value rjg^ and for the second one 
from this r^2 to some r^2 behind the barrier. Both values depend on temperature, of course, as 
may be inferred from Fig.|^. We took the values r^i2 = 0.47, 0.46, 0.45, 0.44, 0.42 and = 
1.05, 0.96, 0.87, 0.74, 0.60 for T = 1, 2, 3, 4, 5 MeV, respectively. The results are put together 
in the following Table: 



T 




hWm 




Vm 


-Cs 


B 






Vs 


MeV 


MeV 


MeV 


MeV/h 




MeV 


MeV 


MeV 


MeV/h 




1 


538 


0.82 


1.8 


1.1 


109 


7.20 


0.50 


1.2 


1.2 


2 


509 


0.94 


4.2 


2.2 


131 


5.95 


0.59 


3.2 


2.7 


3 


472 


0.93 


4.8 


2.5 


134 


4.31 


0.59 


3.8 


3.2 


4 


407 


0.87 


4.6 


2.7 


152 


2.59 


0.60 


4.0 


3.3 


5 


335 


0.79 


4.6 


2.9 


200 


1.18 


0.67 


4.1 


3.0 



Several inferences can be drawn immediately: 

- In the values of rj and (3, as well as in their variation with T there is not much difference between 
the barrier and the minimum. 

- The basic condition on the validity of (p.4|) is satisfied for temperatures up to 3 MeV at best. 

- It is interesting to note that increases with T, albeit only slightly. 

A closer inspection shows the quantity 27rT/{hws) to be definitely larger than unity. This implies 
that collective quantum effects do not contribute to fission dynamics for the examples chosen here. 



According to ||29| for 2nT/{htUs) <^ 1 such effects would change (|3.4| ) to R = fgRx, but quite 



apparently the quantum correction factor fq plays a role only at temperatures smaller than 1 MeV. 
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4 Summary, conclusion and outlook 



In the previous sections we have presented a detailed microscopic study of transport coefficients for 
fission. We concentrated on average motion as in the range of temperatures considered the Einstein 
relation holds true. The latter allows one to deduce the essential diffusion coefficient from friction: 
D = 'jT. Those for average motion, M, 7 and C, have been deduced within the quasi-static picture 
in which linear response theory is applied to describe motion around a local thermal equilibrium. 
Following earlier work, the contribution from the "heat pole" to the response functions has been 
neglected. As described in P], such a restriction is closely related to what in the context of nuclear 
physics one commonly associates with the "adiabatic" picture, in contrast to the "diabatic" one 
one would expect to apply for "fast" motion. 

The response functions have been calculated within a realistic two-center shell model. Effects of 
collisions have been accounted for by using self-energies having both real and imaginary parts. The 
latter are allowed to depend not only on temperature but on frequency as well, in which way memory 
effects of the collision term are simulated. Locally collective motion is treated self-consistently, in 
the sense that the structure of the associated response resembles the one known from the RPA of 
undamped motion. 

Various procedures have been described to deduce transport coefficients from the microscopically 
computed response, which generalize previous descriptions. They have been evaluated as function 
both of temperature as well as of the variation of the shape along the fission path. Considering 
all the effects taken into account, such a study has never been reported before. The results found 
invite us to draw the following conclusions. 

• Perhaps one of the most striking features is the weak dependence of both (3 as well as of rj on the 
shape parameter, as encountered here along the fission path. 

• The temperature dependence of these two quantities is similar to the one reported previously 
in IP] for vibrations of ^^^Pb: they increase with T, to eventually reach some saturation around 
T ^ 4 MeV. 

• This behavior is in agreement with findings reported in [Q], at least qualitatively. Our damping 
rate is somewhat smaller, and its variation with the shape is weaker than needed in for an the 
analysis of the 7-ray multiplicity encountered for fission of '^'^'^Th. 

• However, it must be said that we have not yet striven for quantitative agreements. For instance, 
there is room for a more appropriate choice of the two parameters of coUisional damping. In the 
present paper only their "standard choice" has been used, but as discussed in these values 
are open for changes within a certain margin, and as demonstrated both in |^ and in ||^, such 
modifications will change somewhat (to less than a factor of 2) the values of the transport coefficients 
as well as their variation with T. Furthermore, one should mention the influence of pairing, which 
still might be important at the smaller T value of 1 MeV, but which we have left out for the sake of 
simplicity. Moreover, the influence of angular momentum ought to be taken into account (see below). 
It should be said, of course, that the numbers extracted from comparisons with experiment are model 
dependent. Moreover, in simplified studies often Kramers' picture is applied, even in cases where 
the barrier is too low to guaranty that the flux across the barrier can be described by Kramers' 
solution of the Fokker-Planck equation. This problem can only be cured by performing genuine 
dynamical studies along the hues reported in |]2D[, but where all ingredients of 
the macroscopic descriptions with Fokker-Planck or Langevin equations come from the same theory. 

• Finally, we like to come back to angular momentum, once more, which influences fission dynamics 
in a manyfold way. First of all, it changes the quasi-static energy through the centrifugal potential. 
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Generally, this reduces the height of the barrier and thus restricts further the range of temperatures 
for which ( p.4| ) is applicable. The centrifugal force will also modify the local stiffness and thus the 
coupling constant k through (p.8| ), which in turn effects the local response in various ways. However, 
even in zero frequency limit rotations will have sizable impact on the transport coefficients. This 



has been demonstrated in |^3[ for the case of friction. The reason for such a behavior is found in the 
fact that in the rotating frame the level structure may become very different from the one without 
rotations. 
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Figure 1: The imaginary part of the collective response function (in arbitrary units) as function of 
a;, calculated at the points A and C for T = 1 and 2 MeV. 



19 



Z( Z2 2 

Figure 2: Lower part: The single particle potential V{p, z) along the z-axis; Vq and V are the 
heights of the potential barrier at the origin, calculated with and without the neck-correction term, 
respectively; Upper part: The associated nuclear shape as specified by an equipotential surface. 



Figure 3: The cold liquid-drop energy surface as functions of Zq and 6 for symmetric fission of ^^^T/i. 
The long- and short-dashed curves show the minimal energy fission path and the configuration 
where two fragments start to separate, respectively. ( There exists no postscript file for this figure; 
the authors will be glad to send a hard copy upon request) 
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Figure 4: The free energy as a function of the shape variable ri2 along the fission path for temper- 
atures of 1,2,3,4, and 5 MeV . 
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Figure 5: The local stiffness C(0) of the zero frequency hmit as a function of T at the points A, C 
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Figure 6: The temperature dependence of the friction coefficient Eq.( |2.24D in the zero-frequency 
hmit computed with the oscillator well. The shape variable Q corresponds to 6. The computations 
are performed for the double magic nuclei ^'^^Yb without spin-orbit coupling (upper) and "^^^Pb with 
spin-orbit coupling (lower) 
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Figure 7: The friction coefficient 7 obtained from the oscillator fit, shown as function of T at the 
points A-E 
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Figure 8: The friction coefficient 7 as a function of for T =1-5 MeV 
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Figure 9: The inertia as a function of T at the point A, calculated by three different methods: (i) 
by the oscillator ffi, (ii) by way of Eg. ( |2.18| ) and (iii) after obeying the sum rule (see text). 
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Figure 10: The inertia M satisfying the energy weighted sum rule as a function of ri2 for T —1-5 
MeV 
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Figure 11: P — 7/M as function of ri2 for T —1-3 MeV; for comparison we plot /3 = ^m.w/Mirrot 
and the empirical values from the analysis of Frobrich et. al. 
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